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ABSTRACT 

Recent observations of long and short gamma-ray bursts have revealed a puzzling X-ray ac- 
tivity that in some cases continues for hours after the burst. It is difficult to reconcile such time 
scales with the viscous time scales that an accretion disk can plausibly provide. Here I discuss 
the accretion activity expected from the material that is launched into eccentric, but gravita- 
tionally bound orbits during a compact binary merger coalescence. From a simple analytical 
model the time scales and accretion luminosities that result from fallback in the aftermath 
of a compact binary merger are derived. For the considered mass range, double neutron star 
binaries are relatively homogeneous in their fallback luminosities. Neutron star black hole 
systems show a larger spread in their fallback behaviour. While the model is too simple to 
make predictions about the detailed time structure of the fallback, it makes reasonable pre- 
dictions about the gross properties of the fallback. About one hour after the coalescence the 
fallback accretion luminosity can still be as large as ~ 10 45 erg/s, a fraction of which will 
be transformed into X-rays. Large-scale amplitude variations in the X-ray luminosities can 
plausibly be caused by gravitational fragmentation, which for the high-eccentricity fallback 
should occur more easily than in an accretion disk. 



1 INTRODUCTION 



Recent observations by Swift and Hete II have confirmed the long- 
held suspicion that short GRBs have a different central e ngine than 
long bursts: they occur at systematically low er redshifts dFox et all 
20050, in galaxies with (Hjorth et al. 2005) and without star for- 
mation iBerg eretalJ 120051 ; iBarthelmv etafl 120051) , the emitted 
gamma-ray energy is s ubstantially lower compared to long bursts 
dSoderberg et al .1120060 and no supernova explosion seems to ac- 
comp any the burst teloom et alj200rj ; lHiorth et alj2005MFox et all 
2005). These observations are naturally explained by the most pop- 
ular central e ngine for short bu rsts, the coalescence of a c ompact bi- 



nary system (Blin nikov etai]|l984lEichler et al.| [l989: Paczynski 
1991), but alternative m odels do exist (e.g. iMacFadven et al.l 
' 20051 : iLevan et alj [20061) . The observations, however, have also 



revealed long-lasting X-ray activity after t he main gamma-ray 
emission phase for about half of the bu rsts teurrows et alj|2005l : 
iNousek et ail 12004 lO'Brien et ail l200fjh. Several possible expla- 
nation s have been suggested (e.g. iKing et alj [20051 : IZhang et al.l 
2006), among them is the fragmentation of the outer parts of an 
accretion disk due gravitational instabilities jPerna et al]2006l) . 
This late-time X-ray flaring is considered a challenge to most mod- 
els of short GRBs. In this letter, I discuss the X-ray signature from 
fallback material in the aftermath of a compact binary merger. A 
fraction of the debris material is launched into eccentric orbits and, 
at some later time, falls back towards the central remnant. For this 
material, the time scales are not set by viscous disk time scales, but 



instead by the distribution of eccentricities within the fallback ma- 
terial which is a result of the torques that acted during the merger 
process. This naturally leads to time scales that are orders of mag- 
nitudes larger than usual viscous disk time scales. 



2 A SIMPLE ANALYTICAL MODEL 

The question of fallback time scales is addressed by means of a 
simple analytical model whose initial conditions are taken from 
a set of recent, three-dimensional smooth particle hydrodynam- 
ics (SPH) calculations of the merger of both double neutron stars 
(DNS) and neutron star black hole (NSBH) systems. The details of 
the inp ut physics and the numerical techniques are described else- 
where (Rosswog & Davies 2002; Rosswog & Liebendorfer 2003; 



Rosswog et al. 



2006: 



Rosswog 



20031: iPrice & Rosswod 120061 : iRosswog & Pried 



2006). Here we only summarise the main physi- 



1 It has recently been suggested (Berger et al. 2005) that at least 1/4 of 
the short bursts could lie at redshifts z > 0.7 and would therefore implie 
sunstantially larger isotropised energies than was inferred for the first set of 
detected short GRBs with afterglows. 



cal parameters of the calculations, see Tab. [TJ A few tens of mil- 
liseconds after such a binary encounter the remnant consists of the 
following parts: a central object (either a neutron star-like object 
or a black hole), a disk, some material flung into eccentric orbits, 
but gravitationally bound to the remnant ("fallback material") and 
some dynamically ejected debris material. The details of the geom- 
etry depend on the type of system (DNS or NSBH) and on the mass 
ratio. As an example, Fig. [TJ shows the matter distribution 22.3 ms 
after the merger of a 1.1 and 1.6 Mq DNS system. 
Restricted by the Courant-Friedrichs-Lewy time step criterion (typ- 
ical time steps are of the order 10~ 7 s) the evolution of the merger 
remnant can never be calculated hydrodynamically up to the in- 
teresting time scales of hours. We therefore adopt the following 
simple analytical model. The SPH-particles that are launched into 
eccentric orbits, but are still bound to the remnant, are treated as 
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Table 1. Runs that are used as an input for the analytical model. DNS refers 
to double neutron star systems, NSBH to neutron star black hole binaries. 
Mi, M2 are the masses of the binary components (solar units), q is the 
mass ratio, # part, is the total number of SPH particles in the simulation, 
mfb _2 is the amount of fallback mass in units of 10 -2 Mq and Eg, the 
corresponding energy. 



run 


type 


Mi, M2, q 


#part. 


"Sfb,-2 


l°g(£fb[e 


DA 


DNS 


1.4, 1.4, 1.000 


2 007 516 


3.61 


50.9 


DB 


DNS 


1.1, 1.6, 0.733 


400 904 


3.21 


51.0 


DC 


DNS 


1.1, 1.1, 1.000 


211 296 


3.22 


50.8 


DD 


DNS 


1.2, 1.2, 1.000 


211 296 


2.84 


50.7 


DE 


DNS 


1.3, 1.3, 1.000 


211 296 


3.47 


50.8 


DF 


DNS 


1.5, 1.5, 1.000 


211 296 


3.36 


50.9 


DG 


DNS 


1.6, 1.6, 1.000 


211 296 


3.61 


50.9 


NA 


NSBH 


1.4,4.0, 0.350 


600 581 


7.40 


51.9 


NB 


NSBH 


1.4, 6.0, 0.233 


208 165 


4.57 


51.8 


NC 


NSBH 


1.4, 14.0, 0.100 


2 971 627 


4.05 


51.7 


ND 


NSBH 


1.4, 16.0, 0.088 


1 005 401 


2.91 


51.6 


NE 


NSBH 


1.4, 18.0, 0.078 


1 497 453 


0.015 


49.6 




x [km] 

Figure 1. Contours of the mass density (10 7 — 10 14 g cm -3 ) 22.3 ms after 
the merger of a double neutron star system with 1.1 and 1.6 Mq. In this 
case, about 0.03 Mq are launched into eccentric orbits and will finally fall 
back towards the centre. 



test masses in the gravitational field of the enclosed mass, M, i.e. it 
is assumed that hydrodynamic pressure forces can be neglected and 
the dynamics of each particle can be approximated as a point-mass 
two-body problem. For each of these particles the angular momen- 
tum, Ji, and their energy, Ei, can be determined. Together with 
the particle mass m% <C M these two quantities yield the orbital 
eccentricity 



1 + 



2E t J 2 
G 2 m\M 2 ' 

The semi-major axis 

GMrrii 

di = 

2Ei 



(1) 



(2) 



then provides the maximum and the minimum distance of the par- 
ticle orbit from the origin 



a,i(l 



(3) 



The time until a particle reaches a given radius R^s can be calcu- 
lated by integrating the radial equation 




V(r) 



J? 



2rair 2 



(4) 



where V is the potential energy and the positive (negative) sign 
refers to the motion away from (towards) the origin. The time that 
elapses while a particle moves from a radius r\ to a radius r2 can 
be written as 



r dr 



VAr 2 + Br + C' 



(5) 



where the constants A = B = 2GM and C = - \ have 

been introduced. The integral appearing in eq. {5j can be solved 
analytically 



VAr 2 + Br + C 



( 2Ar + B 



A 




Figure 2. A sample of fallback trajectories resulting from a double neutron 
star merger with 1.1 and 1.6 Mq. The dissipation radius, Rdis m this case 
was 100 km. 



where D = 4AC — B 2 . The time that elapses from a particle's 
present radius, Vi, until its energy is dissipated at a radius i?dis is 
then given by the contributions with the appropriate signs: 



Ti = 



+ Ir. 



i-Rd 



for 
for 



f- > 
•fi <0 



(7) 



For the radius, where the fallback energy is dissipated, 7?dis, 
we take in the double neutron star case the disk radius at the end 
of the numerical simulation. This is conservative in the sense that 
the disk will shrink, see below, and therefore the assumption of a 
radius fixed at i?di s yields shorter timescales and lower accretion 
^ luminosities. For the black hole cases, I choose R^is = WGM/c 2 , 
so just outside the innermost stable circular orbit of a non-rotating 
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DNS 2x1.1 
DNS 1.1 & 1.6 
DNS 2 x 1.2 
DNS 2 x 1.3 
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Figure 3. Fallback accretion luminosity for various compact binaries. Circles refer to double neutron star systems (DNS), triangles to neutron star black hole 
mergers (NSBH). Note that the rightmost point is determined by the mass resolution of the simulation. For reference a straight line with slope 5/3 is shown. 



(Schwarzschild-) black hole at -Risco = 6GM/c 2 . The details of 
the fallback times and energies change slightly with i?dis, but none 
of the conclusions of this paper depends on the exact numerical 
value of i?dis- 

As an illustration, I show in Fig. [2] a set of fallback trajecto- 
ries. The initial conditions were taken from the double neutron star 
merger calculation with 1 . 1 and 1 .6 Mg that is shown in Fig.[2] 200 
randomly chosen trajectories out of the more than 5800 fallback 
particles are plotted. There is a broad distribution of eccentricities 
and fallback times, while some particles return from their initial 
position directly towards the central object, others follow highly 
eccentricity orbits out to distances of many 10 000 km before they 
come to a halt and fall back towards the centre. 



3 FALLBACK ACCRETION 

Fig. [3] shows the accretion luminosities, L acc = dEn>/dt, derived 
for the various DNS and NSBH systems. Here, denotes the 
difference between the potential plus kinetic energy at the start ra- 
dius, n, and the potential energy at the dissipation radius, -Rdis- 
The curves have been obtained by binning the energies contained 
in the fallback material, En,, according to the corresponding fall- 
back times, Ti, see eq. Q. A fraction ex of this energy is channelled 



into X-rays, L x = exiacc 

The double neutron star cases form a rather homogeneous class 
with respect to their fallback accretion, in all cases the fallback 
material is approximately 0.03 Mq, see Tab. Q] After an initial, 
short-lived plateau, the l uminosity falls off with tim e close to the 
expected 5/3-power law jReesll 19881 ; |Phinney|[l98^) . It has to be 
pointed out that the last point in these curves is determined by the 
numerical mass resolution in the hydrodynamics simulations and 
should therefore be interpreted with some caution. All other points 
should be a fair representation of the overall fallback activity. Typ- 
ically, the X-ray luminosity about one hour after the coalescence 
is Lx ~ (fn") ' 10 44 er g/ s - F° r tne investigated mass range, the 
spread in the luminosities one hour after the coalescence is about 
one order of magnitude. 

The neutron star black hole cases show a larger diversity. The mass 
in the fallback material of different mass ratios varies by about a 
factor of 500, see Tab. [T] an hour after the merger the accretion 
luminosities of the different NSBH systems differ by about two or- 
ders of magnitude. Also the involved time scales change strongly 
with the binary mass ratio. For example, the 1.4 and 4 Mq NSBH 
case does not produce much eccentric fallback material. Accretion, 
at least to the resolvable level, is over in ~ 0.2 s. This accretion 
period may produce a short GRB, but probably not much X-ray ac- 
tivity. The 1.4 and 18 Mq NSBH system is at the other extreme: its 
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peak luminosity is lower by three orders orders of magnitude but 
extends (at a resolvable level) up to about one hour. The mass ra- 
tios inbetween could possibly produce a (weak) GRB and extended 
X-ray activity up to about one day after the burst. 



4 SUMMARY 

A simple analytical model was presented and used to derive the 
fallback time scales and accretion luminosities in the aftermath of 
a compact binary merger. The initial conditions are taken from a 
large set of 3D hydrodynamical simulations. The main assumption 
of this model is that the fallback material is ballistic, i.e. that hydro- 
dynamical pressure forces can be neglected. This is a good approx- 
imation for the largest part of the time spent along the individual 
trajectories, therefore, the fallback time scales can be considered 
valid estimates. The assumption may only break down when the 
vicinity of remnant is approached, but this will only cause minor 
corrections to the fallback time scale and luminosity estimates. 
For the range of mass ratios (1.1 to 1.6 M©) that I have considered, 
double neutron star systems show a rather homogeneous fallback 
behaviour. At one hour after the coalescence their X-ray luminos- 
ity is still ~ 10 44 (ex/0.1) erg/s. Neutron star black hole mergers 
provide a larger spread in their fallback behaviour. An hour after 
the tidal disruption of the neutron star, their range in fallback lumi- 
nosities is about two orders of magnitudes. Some mass ratios could 
produce a GRB, but not much X-ray activity, others may produce 
X-ray flaring only but no GRB, and yet other systems can produce 
both. 

Since this simple approach does not account for self-gravity and 
gas pressures, it tends to produce too smooth results. The shown 
luminosity curves should therefore be interpreted as temporal av- 
erages. The original accretion disks, such as the one seen in Fig.Q] 
extending out to ~ 200 km, are consumed on a viscous time scale 
of 

where a is the Sha kura-Sunyaev viscosity parameter 
dShakura & Sunvaevl Il973l) . R the disk radius, Qk the Ke- 
pler frequency and M c . . is the mass of the central object. The 
outer parts of this disk are constantly fed by material that continu- 
ously falls back so that a low mass remnant disk is maintained until 
late times. The observed X-ray activity is expected to be produced 
by this disk and its interaction with the high-eccentricity fallback 
material. 

Large-scale amplitude variations in the lightcurves ("flares") can 
occur if the fallback material becomes subject to self-gravitational 
instabilities. Whether a lump of gas is prone to this instability 
depends on the ratio o f the coo l ing tim e, r coo i, and the local 
dynamical time, Td yn . iGammi i l l200ll) finds for the case of 
thin accretion disks the condition r coo i ^ 3rd yll . Gravitational 
fragmentation i s common in the ou ter regions of the accretion 
disks of AGNs (Shlosman et al. 1990) and of young stellar objects 
dAdams&Lmlll993l) . iPerna et alj ( 120061) have suggested that such 
a fragmentation in the outskirts of an accretion disk could be the 
common mechanism behind the flaring observed in both long and 
short bursts. They note, however, that it may be difficult to reach 
the suitable conditions in the case of a compact binary merger. 
The frag mentation idea has r ecently been examined in more 
detail by IPiro and Pfahl (2006) for the case of collapsar disks. 
They argue that fragmentation may occur near the radius where 



alpha-particles are photodissociated since this consumes about 7 
MeV of thermal energy per nucleon and thus provides an efficient 
cooling mechanism. They suggest that even neutron star-like 
objects of low masses may form, but conclude that their disruption 
would be too rapid to account for the observed, post-GRB X-ray 
flaring. 

Some the eccentric fallback material could undergo grav- 
itational fragmentation into one or several bound objects. 
Since there is a large spread in orbital eccentricities, it seems 
unlikely that a large fraction of the fallback material ag- 
glomerates into a single fragment. If fragmentation occurs, 
it will do so only locally and for a small fraction of the ma- 
terial. A fragment of mass nif will release a flare energy of 

*~ lO^ergfe) (f^) (g^) . When such 

a fragment of radius Rf falls back towards the centre, it becomes 
tidally disrupted at a distance i? t id = (Mco./rrif) 1 / 3 ■ Rf, or, 
if the remnant disk radius is larger than Rtid, the fragment will 
impact with nearly free-fall velocity on that disk. In both cases, it 
will trigger violent X-ray flaring activity. 

In this picture the GRB is produced by the massive accretion 
disk that can form early on in the merger. The GRB duration 
is determined by the viscous time scale on which this disk is 
consumed. The late-time X-ray activity, in contrast, is caused by 
the fallback of a smaller amount, about a tenth of the disk mass, 
of high-eccentricity material which is a natural outcome of the 
merger. The fallback time scales are set by the distribution of 
eccentricities in this material, which is a result of the torques that 
acted during the merger process. 



ACKNOWLEDGEMENTS 

The calculations in this paper have partly been performed on the 
JUMP computer of Forschungszentrum Jiilich. 



REFERENCES 

Adams F. C, Lin D. N. C, 1993, in Protostars and Planets in, 

Levy E. H., Lunine J. I., eds., pp. 721-748 
Barthelmy S. D., et al., 2005, Nature, 438, 994 
Berger E. et al., 2005, Nature, 438, 988 
Berger E. et al., 2006, submitted to Ap J, |astro-ph/06 1 1 1 28 1 
Blinnikov S. I., et al., 1984, Soviet Astron. Lett., 10, L77 
Bloom J. S. et al., 2006, ApJ, 638, 354 
Burrows D. N. et al., 2005, Science, 309, 1833 
Eichler D., Livio M., Piran T, Schramm D. N., 1989, Nature, 340, 

126 

Fox D. B. et al., 2005, Nature, 437, 845 
Gammie C. F, 2001, ApJ, 553, 174 
Hjorth J. et al, 2005, Nature, 437, 859 
Levan, A.J. et al . , j astro-ph/060 1 3 32 1 (2006) 
MacFadyen A. I., Ramirez-Ruiz E., Zhang W., 2005, 
|astro-ph/0510192| 

Nousek J. A. et al., 2006, ApJ, 642, 389 

O'Brien P. T. et al., 2006, ApJ, 647, 1213 

Paczyfiski B., 1991, Acta Astron., 41, 257 

Perna R., Armitage P. J., Zhang B., 2006, ApJL, 636, L29 

Piro A., Pfahl, E., 2006, submitted to ApJL, |astro-ph/06 1 0696 1 

Phinney E. S., 1989, in IAU Symp. 136: The Center of the Galaxy, 

Morris M., ed., pp. 543 
Price D., Rosswog S., 2006, Science, 312, 719 



Fallback accretion 5 



King, A.R. et al., 2005, ApJ 630, LI 13 
Rees M. J., 1988, Nature, 333, 523 

Rosswog S., Revista Mexicana de Astronoma y Astrofsica, 2006, 
in press 

Rosswog S., Davies M. B., 2002, MNRAS, 334, 481 
Rosswog S., Liebendorfer M., 2003, MNRAS, 342, 673 
Rosswog S., Price D., 2006, in prep. 

Rosswog S., Ramirez-Ruiz E., Davies M. B., 2003, MNRAS, 345, 
1077 

Shakura N. I., Sunyaev R. A., 1973, A & A, 24, 337 
Shlosman I., Begelman M. C, Frank J., 1990, Nature, 345, 679 
Soderberg A. M. et al., 2006, ApJ, 650, 261 
Zhang B. et al., 2006, ApJ, 642, 354 



